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SUMMARY 


Heat-transfer rates have been measured in free flight along the 
stagnation line of an unswept cylinder mounted transversely on an axial 
cylinder so that the shock wave from the hemispherical nose of the axial 
cylinder intersected the bow shock of the unswept transverse cylinder. 
Data were obtained at Mach numbers from 2.53 to 5*50 and at Reynolds 
numbers based on the transverse cylinder diameter from 1.00 x 10^ to 
I.87 x 10 6 . 

Shadowgraph pictures made in a wind tunnel showed that the flow 
field was influenced by boundary- layer separation on the axial cylinder 
and by end effects on the transverse cylinder as well as by the inter- 
secting shocks. Under these conditions, the measured heat-transfer rates 
had inconsistent variations both in magnitude and distribution which 
precluded separating the effects of these disturbances. The general 
magnitude of the measured heating rates at Mach numbers up to 3 was from 
0.1 to 0.5 of the theoretical laminar heating rates along the stagnation 
line for an infinite unswept cylinder in undisturbed flow. At Mach 
numbers above A the measured heating rates were from 1. 5 to 2 times the 
theoretical rates. 


INTRODUCTION 


Along with the increasing speed of high-performance aircraft has 
come, necessarily, an increasing emphasis on the problem of aerodynamic 
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heating. In recent years research has been devoted to determining the 
heat transfer to blunt bodies, such as cylinders, and wing leading edges 
subjected to uniform high-speed flow in order to obtain information use- 
ful in design. The problem of determining the heat transfer to wings or 
bodies subjected to a nonuniform flow field, as in the case where a wing 
or body is subjected to shocks and other disturbances originating from 
another body, has more recently come under scrutiny with the advent of 
such vehicles as the X-15 and Dyna-Soar. 

To investigate the influence of a nonuniform flow field on heating 
rate to an unswept cylinder, a two-stage rocket model system was designed 
and launched from the NASA Wallops Station. The cruciform- shaped test 
configuration, which was mounted on the nose of the second stage, con- 
sisted of two cylinders having nearly equal diameters and instrumented 
to measure skin temperatures along the stagnation line of the transverse 
cylinder. The shock wave from the hemispherical nose of the axial cyl- 
inder impinged on the bow shock of the unswept transverse cylinder. Also, 
the upstream axial cylinder and its boundary layer caused other disturb- 
ances in the flow on the transverse cylinder. Aerodynamic heating rates 
were obtained along the stagnation line of the unswept transverse cylin- 
der for Mach numbers up to 5. 50 and Reynolds numbers based on the trans- 
verse cylinder diameter from 1.00 X 10^ to I .87 X 10^. Results from 
this test and some concluding remarks are presented herein. 


SYMBOLS 


D diameter of transverse cylinder, 0 . 0625 ft 

d diameter of hemispherical-tipped axial cylinder, 0.0754 ft 

h measured heat-transfer coefficient, Btu/sec-ft^-°R 

h^ theoretical heat-transfer coefficient for laminar flow, 

Btu/ft 2 -sec-°R 

k coefficient of thermal conductivity, Btu-ft/sec-ft^-°R 

l length along stagnation line between thermocouple 

stations, ft 

M Mach number 

p pressure, lb/sq ft 

Q spanwise conduction at station it, Btu/sec-ft^ 
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q heating rate, Btu/sec-ft 2 

Boo free-stream Reynolds number per foot 

Rp free-stream Reynolds number based on transverse cylinder 

diameter 

r r radius of transverse cylinder, ft 

T temperature, °R 

TC abbreviation for thermocouple in figures 3 and 10 

t time, sec 

V velocity, ft/sec 

x distance along transverse cylinder from model center 

line, ft 

y distance along axial cylinder from hemispherical nose 

tip, ft 

0 circumferential angle, radian 

p density, slugs/cu ft 

t skin thickness, ft 

Subscripts: 

1 inner surface 

n nth term in numerical series (l, 2, 3> • • • , n) 

o outer surface 

t stagnation 

w vail 


00 


free stream 
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MODEL, INSTRUMENTATION, AND TEST 


Model 

The model consisted of a flare-stabilized JATO, 1.52-KS-33>550> 

SM19 (Recruit) rocket motor with the test nose and instrument section 
mounted on the forward end of the motor. The test nose which was a 
cruciform configuration consisted of a hemispherical-tipped cylinder to 
serve as a shock producer, and a transverse cylinder with zero sweep on 
which the heat-transfer measurements were made. The transverse cylinder 
was located so that at high speeds the shock from the hemispherical tip 
would intersect the shock wave of the transverse cylinder. Photographs 
of the model and the test nose are shown in f gures 1 and 2 , respectively. 
The shock producing tip of the test nose and bhe transverse cylinder were 
made of Inconel. The remaining part of the cruciform test nose was made 
of mild carbon steel (SAE 1018). In order to obtain a fast temperature 
response through the wall of the transverse cylinder, a nominal wall 
thickness of 0.032 inch was used. As predicted by use of reference 1, 
this thickness would be sufficient to withstand the heating through the 
peak Mach number of the test. This prediction did not include effects 
due to a nonuniform flow field. 


Instrumentation 

The model was equipped with a telemetering system which transmitted 
the data to a ground receiving station. Three channels were used to 
transmit the accelerations (longitudinal, transverse, and normal) expe- 
rienced by the model, and two channels were used to transmit the temper- 
atures measured along the cylinder semispan an the thermocouple locations. 
Figure 3 shows the seven evenly spaced thermocouples mounted along the 
stagnation line on the inside surface of the transverse cylinder. Tem- 
perature data from thermocouples 1 , 2 , 3 , 5 , o, and 7 were commutated 
and transmitted on one of these temperature channels so that each thermo- 
couple recorded about every 0.1 second. Temperature data from thermo- 
couple 4 were transmitted continuously on the other temperature channel. 
The maximum probable error of the temperature measurement was ±2 percent 
of the calibrated full-scale value, or, in this case, ±40°. 

In addition to the internally carried instrumentation, the model 
was tracked by a CW Doppler velocimeter to provide the velocity-time 
history, and by two different radar sets, a N4SA modified SCR-584 
tracking radar and an AN/FPS-16 radar set. Each provided flight-path 
data for the test. Atmospheric and wind conditions were determined by 
means of a radiosonde launched near the time of flight and tracked by 
a Rawin set AN/GMD-1A. 
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Test 

The model was launched at an elevation angle of 60°. The booster, 
a fin- stabilized M5 JATO rocket motor (Nike), accelerated the model to 
a Mach number of about 3-0. A picture of the model and booster on the 
launcher is shown in figure 4. After a short coast period, the sus- 
tainer rocket ignited, effected separation from the spent booster, and 
accelerated the model to its peak Mach number of 7*4* 

At a Mach number of 5*9^-^ which occurred 1. 14 seconds after ignition 
of the sustainer, the telemeter signal was lost. However, radar tracking 
indicated that the model continued accelerating along its anticipated 
flight path until after sustainer burnout. The telemetered data showed 
that at the time of signal failure the temperatures on the stagnation 
line of the transverse cylinder were rising very rapidly and were 
approaching the structural thermal limit of Inconel. (it may be noted 
that If boundary- layer transition occurred on the transverse cylinder, 
the stagnation line would not necessarily be the hottest location. ) 

At sustainer ignition the model experienced an oscillation in pitch 
which had damped to essentially zero degrees at the time of signal fail- 
ure 1. 14 seconds later. The high temperatures and air loads encountered 
probably caused a breakup of the model nose and loss of telemeter signal 
without completely destroying the vehicle. 

Heating data presented in this report are for the time of flight 
up to 5*34 seconds (Moo = 5 * 50 )> during which time telemeter signals 
were received. This Mach number was reached at an altitude of 10,500 feet 

with a corresponding Reynolds number of I .87 X 10^ based on the trans- 
verse cylinder diameter. 


DATA REDUCTION 


Trajectory Data 

Time histories of model velocity, Mach number, Reynolds number per 
foot, and altitude are presented in figure 5- Variations of the free- 
stream air properties (temperature, density, and pressure) with flight 
time are given In figure 6. The model velocity was obtained by inte- 
gration of the longitudinal acceleration up to the time of telemeter 
failure and thereafter from the CW Doppler velocimeter data. The other 
data of figures 5 6 were obtained from ground radar and radiosonde 

measurements. 
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Heating Rates 

Faired curves showing the variation of skin temperatures with time 
are plotted in figure 7- The curves for different stations end at 
different times because the thermocouples, except for the one at sta- 
tion by were commutated on a single channel. The faired curves of meas- 
ured inside wall temperature were used to compute the outside surface 
temperature. The outside surface temperatures were then used to compute 
one -dimensional heating rates. The computational methods used were those 
of reference 2, which considers a thermally thick wall without curvature 


and with no heat flow parallel to the wall anc: no plate curvature. Effect L 
of wall curvature was not included but would increase the heating rates 8 
less than 5 percent. Faired curves of the experimental one -dimensional 7 
heating rates for two representative stations are shown in figure 8. 9 


The aerodynamic heat-transfer rate at a measurement station is con- 
sidered to be equal to the one-dimensional heating rate plus the rates 
of heat loss at the station due to conduction and radiation. Sample 
computations of radiation showed it to be always less than 2 percent of 
the one -dimensional heating rate and it was therefore disregarded. Like- 
wise, conduction along the skin normal to the stagnation line was dis- 
regarded since estimates showed it to be always less than 5^ percent of 

the one-dimensional heating rates. The estimates were made by use of 
temperature distributions derived from the laminar heat-transfer dis- 
tributions on an unswept cylinder in reference 3. The rate of heat 
accumulation at a station due to spanwise conduction along the stagnation 
line was computed by a finite-difference method by use of the measured 
temperatures and the following equation: 

= k ( T n-l - T n )r(r’d9) k(T n - T n+1 )r(r’de) 
z z ( r ' ae ) z:.(r'ae) 

* ^(Tn-l - T n) - ( T n - T n+l) 

The one-dimensional heating rates were corrected by this spanwise con- 
duction term to determine aerodynamic heating rate. 
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RESULTS AND DISCUSSION 


Shock Position 

An estimation was made of the location of the intersection of the 
shock wave from the hemispherical nose and the bow shock of the trans- 
verse cylinder since it was expected that this location would have a 
strong influence on the distribution of heating along the transverse 
cylinder. In estimating the position of these shocks, use was made of 
unpublished results from tests of a similar nose configuration made in 
the preflight jet at NASA Wallops Station at a Mach number of 2 and in 
the Langley Unitary Plan wind tunnel at Mach numbers of 2.65 an( i 
In those tests the transverse cylinder was about 2 diameters longer in 
span than in the present flight test. Figure 9 presents a shadowgraph 
from each of these tests. In each, a dark line extends from the inter- 
section of the two strong shocks towards the surface of the transverse 
cylinder. Since the flow in this region near the stagnation line of the 
transverse cylinder is subsonic, this dark line does not indicate a 
shock wave but may be a slip plane or the three-dimensional photographic 
effects of the shock intersection. 

Separation of the boundary layer on the axial cylinder in front of 
the transverse cylinder is also apparent in the shadowgraphs. The pres- 
ence of spanwise flow and the general complexity of the flow field around 
the configuration are indicated by the spanwise curvature of the bow 
shock outboard of the shock intersection and the irregularity of the 
shock inboard of the intersection. 

Measurements of the shock shape from these shadowgraphs are shown 
in figure 10 along with shock shapes for several Mach numbers which were 
predicted by the method of reference k. A comparison of the measured 
shock shapes shows very close agreement with the predicted shapes. The 
approximate spanwise location of the intersection of the nose shock wave 
and the bow shock wave is shown in figure 11 as a function of Mach number. 
The location is given in terms of either the nose diameter or the diam- 
eter of the transverse cylinder. Both plots are included since the shape 
of the nose shock (and thus the location of the intersection) is dependent 
on nose diameter; whereas, in subsequent figures spanwise location on the 
transverse cylinder is presented in terms of the diameter of the trans- 
verse cylinder. 


Heat Transfer 

The spanwise distributions of the experimental heating rates on the 
transverse cylinder are presented in figure 12 for several times during 
the flight test over the Mach number range from 2.25 to 5*50. Both the 
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one -dimensional and the aerodynamic heating rites are shown, their dif- 
ference being the correction for spanwise conduction. The outside wall 
temperatures at the thermocouple locations are shown, connected by straight 
lines which indicate the spanwise temperature distributions assumed in 
the computation of the conduction corrections. Also shown on each dis- 
tribution is the spanwise location of the intersection of the nose shock 
and the bow shock of the transverse cylinder as determined from figure 11. 

The theoretical laminar heating-rate distributions were determined 
from the following relation: 

1 - »L( T t - V,) 

with the use of the plotted experimental values of ^wo and theoretical 
values of h^ for the stagnation line of an unyawed infinite cylinder 
(i.e., no variation in h^ with spanwise location) as determined from 

equation (B5) of reference 3* This theory was shown in reference 3 to 
agree within 8 percent with wind-tunnel results at M =■ 4.15 over the 
Reynolds number range of the present test. 

The theoretical turbulent distributions shown in figure 12 were 
determined by use of theory presented in reference 5 and are for a loca- 
tion 5° chordwise from the stagnation line on an unyawed infinite cyl- 
inder at the flight test conditions. They are shown primarily for com- 
parison with the theoretical laminar rates to indicate the predicted 
increase in heating rate for turbulent flow near the stagnation line. 
According to the theory of reference 5j turbulent rates for a given 
chordwise location are a multiple of the laminar rates for that location, 
with the multiplier being a function only of she local Reynolds number 
based on surface length from the stagnation line. Local Reynolds numbers 
at a 5° location for the flight conditions were determined by use of 
Newtonian pressure distribution and the compressible-flow relations of 
reference 6. An error of about 1 percent was introduced by basing the 
theoretical turbulent rates on laminar theory for the stagnation line 
rather than on laminar theory for the 5° location. 

The shadowgraph pictures of figure 9 show that the flow field was 
influenced by boundary-layer separation and by end effects from the 
transverse cylinder as well as by the intersecting shocks. (To deter- 
mine the effect of these intersecting shocks was the original intention 
of the test. ) Figure 12 shows that under these conditions, the measured 
heat-transfer rates had inconsistent variations both in magnitude and 
distribution which precluded separating the effects of these flow dis- 
turbances. For example, no consistent effect of the shock intersection 
on heat transfer was apparent. 
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The relatively high heating rates at the two inboard stations at 
times prior to 3*66 seconds, which caused the temperatures at these two 
stations to rise rapidly (see fig. 7(a)), are apparently not associated 
with Mach number and Reynolds number since they persisted only part way 
through the period of relatively constant flight conditions between 3- 0 
and 4. 5 seconds. The heating rates at the outboard station are, except 
at the two highest Mach numbers, higher than at the neighboring inboard 
stations; this may be due to end effects on the transverse cylinder. 

At Mach numbers up to 3*23 (figs. 12(a) to 12(2)) the experimental 
heating rates were generally much less than laminar theory for an infi- 
nite cylinder, except at the two inboard stations where the previously 
noted high heating rates were sometimes as high as the turbulent theory 
for the 5° chordwise location. At Mach numbers above 3*23 the measured 
rates generally exceeded the laminar theory at all measurement stations. 

For convenient comparison of the experimental data with the laminar 
inf ini te- cylinder theory, the experimental heating rates divided by the 
corresponding laminar theoretical rates, as obtained from figure 12, are 
plotted in figure 13. (Note that the measured heating rate divided by 
the theoretical laminar heating rate is identical to the measured heat- 
transfer coefficient divided by the theoretical laminar heat-transfer 
coefficient since the experimental values of T t - T wo were used to 

compute the theoretical laminar heating rate. ) Except at the two inboard 
stations, the general level of h/h^ along the transverse cylinder was 

between 0.1 and 0. 5 at Mach numbers up to 3- 23. As the Mach number 
increased from 3-23 to 5-50, the general level of h/h^ increased to 

values between 1. 5 and 2.0, with considerable spanwise variation of the 
ratio in most cases. 

At times from 3.84 seconds to 5 *16 seconds the magnitude of h/h^ 

is high near the shock-intersection location relative to values at 
neighboring locations. However, this condition is not evident after 
5.16 seconds, nor before 3*8^ seconds, even though Mach number and 
Reynolds number were approximately constant from about 3*0 seconds to 
4.5 seconds. 

Results of the present free-flight test show a considerably lower 
level of heat transfer along the stagnation line of an unswept cylinder 
in a flow field disturbed by an upstream body than did data obtained 
recently in the Langley Unitary Plan wind tunnel. (See ref. 7* ) Also, 
the present flight test did not show an increase in heat transfer near 
the shock intersection as did the wind-tunnel tests. There were some 
differences between the flight test and the wind-tunnel tests, both in 
geometry and test conditions, which precluded a good comparison of the 
results. 
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CONCLUDING REMARKS 


Heat-transfer rates have been measured in free flight along the 
stagnation line of an unswept cylinder mounted transversely on an axial 
cylinder so that the shock wave from the hemispherical nose of the axial 
cylinder intersected the bow shock of the unswept transverse cylinder. 
Data were obtained at Mach numbers from 2.53 to 5*50 and Reynolds numbers 
based on the transverse-cylinder diameter from 1.00 X 10^ to I .87 X 10^. 

Shadowgraph pictures made in a wind tunneL showed that the flow 
field was influenced by boundary -layer separation on the axial cylinder 
and by end effects on the transverse cylinder as well as by the inter- 
secting shocks. Under these conditions, the measured heat-transfer 
rates had inconsistent variations both in magnitude and distribution 
which precluded separating the effects of these disturbances. The gen- 
eral magnitude of the measured heating rates at Mach numbers up to 3 was 
from 0.1 to 0.5 of theoretical laminar heating rates along the stagnation 
line for an infinite unswept cylinder in undisturbed flow. At Mach num- 
bers above k the general magnitude of the measured heating rates was 
from 1.5 to 2 times the theoretical rates. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , June 22, 1961 . 
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(b) Stations 5; 6, and 7* 
Figure 7 . - Concluded. 
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(a) Station 2. 

Figure 8. - Typical variation of one -dimensional aerodynamic heating 

with flight time. 
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(b) Station 7- 
Figure 8.- Concluded* 
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Figure 9 .- Shadowgraphs of typical shock patterns produced by a similar test configuration. 






Figure 10. - Shock-wave geometry. 
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Figure 12.- Continued. 
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Laminar theory (re'.3) 

Turbulent theory (ref-6) 

O One -dimen sional q 

• Aerodynamic q 

A Wall temperature 

* Approximate location of 

shock intersection 




(h) t = 3 - 84 ; = 2.98; r d = 2.19 x 10 6 . 
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(i) t = 4.03; = 2.97; r d = 1.15 X 10 6 . 


Figure 12. - Continued 
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